Desynchronized wave patterns in synchronized chaotic regions of coupled map lattices 
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We analyze the size limits of coupled map lattices with diffusive coupling at the crossover of low- 
dimensional to high-dimensional chaos. We investigate the existence of standing-wave-type periodic 
patterns, within the low-dimensional limit, in addition to the stable synchronous chaotic states 
depending upon the initial conditions. Further, we bring out a controlling mechanism to explain 
the emergence of standing-wave patterns in the coupled map lattices. Finally, we give an analytic 
expression in terms of the unstable periodic orbits of the isolated map to represent the standing-wave 
patterns. 
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I. INTRODUCTION 

Coupled dynamical systems often arise in nature 
whenever a collective or cooperative phenomenon is 
favoured 0,0)0,0]. In particular, the coupled map lat- 
tice with diffusive coupling (CML) provides a prototype 
model to study various features associated with the co- 
operative evolution of constituent systems 0, 0, 0, 0] • 
One of the important properties of such CMLs is that 
they exhibit size instability, that is, there is a criti- 
cal size on the number of constituents for which stable 
synchronous chaotic state exists. Increasing the num- 
ber of constituents beyond this limit leads to the oc- 
currence of spatially incoherent behaviour (eg., high- 
dimensional chaos). For example, Bohr and Chris- 
tensen [To| have studied such size instability behaviour in 
a two-dimensional coupled logistic lattice. Similar desyn- 
chronization has been found in arrays of couple d sy stems 
represented by nonlinear oscillators 0, ITTl If3 . Il3| . The 
stability of synchronous chaos in coupled dynamical sys- 
tems plays an important role in the study of pattern for- 
mation, spatiotemporal chaos, etc. pi IToL [Til IbH Il5| . 

In general, these studies on size instability are valid in 
most situations. However, we have noted that in certain 
circumstances there is an ambiguity in dealing with these 
systems below the critical system sizes. To be specific, 
there exist certain nontrivial ranges of initial conditions 
for which the CML admits spatial and temporally peri- 
odic solutions in contrast to the usually expected stable 
synchronous chaos. In this brief report, we show numer- 
ically the coexistence of such periodic states with the 
stable synchronous chaotic state well below the critical 
system size and explain the underlying mechanism. 



II. SIZE INSTABILITY IN COUPLED MAP 
LATTICES WITH DIFFUSIVE COUPLING 

Consider an one-dimensional coupled map lattice with 
nearest- neighbour diffusive coupling 0, 0] 

<+i = /(4) + * [/(4 _1 ) + /(4 +1 ) - 2/(4)] . (!) 



where j (= 0, 1, 2, • ■ • , L — 1) represents the lattice sites 
and L is the system size, subject to periodic boundary 
conditions. 

The stability analysis of the synchronized chaotic state 



defined by s n 



„L-l 



in the above CML 



using the procedure derived originally for coupled oscil- 
lators by Heagy et al gives the relation connecting 
transverse Lyapunov exponents (TLEs), in terms of the 
Lyapunov exponent of single (isolated) map, A as 



A fc = A + In 



, , 7rfc 
1 - 4esm I — 



(2) 



The synchronous state is stable only if the TLEs (A fe , 
k = 1, 2, ■ • • L — 1) are all negative. 

The above relation J5J) can also be obtained by means 
of a direct perturbation of the form 



c n = s n + S exp {^—j-^j ex P i^ 71 ) ' 



(3) 



as considered by Bohr and Christensen 0] ■ 

The synchronous state loses its stability when the long 
waves (lowest mode) are unstable [13 . This means that 
for A 1 > the synchronous state is unstable. Thus, 
substituting A 1 = in Eq. (0, one obtains the maxi- 
mum/critical lattice size (L c ) that supports stable syn- 
chronous state as 



L r = hit 



sin 



1 / / l-e- A o 



(4) 



Now, let us consider a coupled logistic lattice with diffu- 
sive coupling (CLL) where each lattice site in Eq. (JIJ is 
occupied by the logistic map 

f(x)=tix(l-x), 16(0,1), Me (0,4). (5) 

In particular for the choice fi = 3.5732, the Lyapunov 
exponent of single (isolated) map is positive (i.e., A ~ 
0.057 > 0) and chaotic. In this case, for coupling strength 



2 



e = 0.2, the critical lattice size (L c ) is found to be 11 from 
Eq. That is, upto the lattice size 11 the CLL exhibits 
synchronous chaos and for lattice size > 12 the synchro- 
nization is found to be lost, thereby confirming the size 
instability in the diffusively coupled logistic lattices. 



III. EXISTENCE OF MULTIPLE STABLE 
STATES IN COUPLED MAP LATTICES 

When random initial conditions are assumed, in most 
cases, the above CLL with diffusive coupling exhibits 
stable synchronous chaos for 2 < L < 11 as predicted 
by Eq. However, there are certain ranges of initial 
conditions for which CLL shows some interesting asyn- 
chronous spatiotemporal patterns even for L < L c . For 
example, for L = 8 and with the choice of initial con- 
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FIG. 1: The possible spatiotemporal patterns in coupled lo- 
gistic lattice and JjjJ with L = 8 for different initial con- 
ditions, (a) single standing wave, (b) double standing waves 
(c) single standing wave with a temporal period-2 orbit in one 
half and period-4 orbit in the other half of the lattice and (d) 
synchronized chaos. 

ditions {x J }f^ = (O- 1 ' °- 01 > °' 7 ' °- 2 > °' 65 > O- 1 ' °- 15 > 
0.001} or several nearby initial conditions, the CLL ex- 
hibits a standing wave type pattern as shown in Fig.^a). 
But for the same lattice size, if we choose a different set 
of initial conditions, {x J }f^ = {0.0004, 0.0001, 0.0003, 
0.0005, 0.00045, 0.00018, 0.00016, 0.00001} or the nearby 
points, two standing waves with different amplitudes are 
produced within the CLL as in Fig. GJb). Similarly a 
disturbed standing wave pattern as shown in Fig. ^c) 
is possible to be exhibited by the CLL for many choices 
of initial conditions. A synchronous chaos, as one would 
expect for L = 8 from the theory, is also exhibited by 
the CLL as depicted in Fig. ^d), for most of the ran- 
dom choices of initial conditions. This kind of multiple 
stable solutions is also observed in the CLL for other 
lattice sizes, namely L = 6, 7, 9 and 10 which are also 
less than L c . The occurrence of various spatiotemporal 



patterns for different lattice sizes of the CLL and their 
percentage of occurrence are shown in Table [I] In order 
to quantify the percentage of occurrence of different spa- 
tiotemporal patterns, we have used 10 6 sets of random 
initial conditions (i.c.'s) in the interval to 1 and identi- 
fied the number of i.c.'s which lead to a specific pattern 
as indicated in Table [fl We have further confirmed our 
assertion by analyzing the same systems in different com- 
puting environments such as Intel Pentium 4, Sun Sparc 
server/workstation and Compaq Alpha workstation. 

In addition, by considering the whole CML of size L 
as a single L-dimensional map we have verified that the 
above spatiotemporal periodic structures are essentially 
the stable fixed points of this map. For example, the pe- 
riodic structure in Fig.^a) represents a fixed of point of 
period-2 of the eight-dimensional map and the eigenval- 
ues of the corresponding Jacobian matrix all are having 
magnitude less than unity. In a similar fashion one can 
verify that all the periodic structures are the stable fixed 
points of corresponding periods. Thus, in addition to the 
stable synchronized manifold, there exists other invariant 
sets corresponding to stable periodic structures. 

TABLE I: The possible spatiotemporal patterns for different 
lattice sizes and their percentage of occurrence for the cou- 
pled logistic lattice sampled over a set of 10 6 random initial 
conditions (IC's). The parameters are fixed as /i = 3.5732 in 
Eq. and e = 0.2 in Eq. Q. Here, SW: standing wave, 
sync: synchronized. 
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IV. EMERGENCE OF STANDING WAVE 
PATTERNS BY CONTROLLING 

The extraordinary behaviour of the CML with diffu- 
sive coupling showing standing wave patterns well below 
the critical lattice size (L c ) can be explained as follows. 
The second term in the right hand side of Eq. can be 
considered as a kind of force or perturbation applied to 
every lattice point in the CML and we call it as the cou- 
pling force. In fact this force on a particular lattice point 
is developed either due to a mismatch in the parameters 
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of the neighbouring lattice points or due to differences in 
their initial conditions or both. If the neighbouring lat- 
tice points are identical then this force is formed due to 
variation in the initial conditions and our system indeed 
falls under this category. In general, the strength of the 
coupling force at all the lattice points approaches zero 
when they oscillate towards synchronization with their 
neighbours, and usually this will happen for L < L c . 
But in the case of L > L c , this force at every lattice 
point oscillates periodically or in a chaotic manner, giving 
rise to various spatiotemporal patterns, including stand- 
ing waves. However, as we have pointed out above that 
under certain circumstances, (i.e., for certain ranges of 
initial conditions), even for L < L c , the coupling force of 
each and every lattice point oscillates periodically with 
different amplitudes and thereby makes the CML to ex- 
hibit spatiotemporal periodic (standing wave) solutions. 
In general, the periodic oscillations in the coupling force 
may be of any period and this fixes the number of stand- 
ing waves produced within the lattice. 

In order to understand the mechanism behind the 
emergence of spatiotemporal periodic structure, let us 
now consider a specific case of the periodically oscillating 
coupling force of period two, that is, the force oscillating 
between two fixed amplitudes, say, k\ and ki so that a 
single standing wave is formed in the CML. Then the am- 
plitudes of the coupling force at the j th lattice site will 
alternate between the numerical values k\ and k? 2 , where 
j = 1, 2, . . . , L. Thus, the evolution of j th map in the 
CML can be effectively described by the equation 

x n+l = f( x n) + x ii+2 = f( X n+l) + ^2 • (6) 

For a given lattice point j, this is nothing but a single 
logistic map with a periodic kick of period-2 (modified 
map). It is now obvious to note from Eqs. I© and 
that the original coupled map lattice exhibiting a stand- 
ing wave pattern can be decomposed into L number of 
modified maps such that the dynamics of Q is essentially 
mimicked by the set JSJ). Thus studying the evolution of 
L decoupled modified maps I© with allowed sets of val- 
ues for k\ and k J 2 is equivalent to that of the original 
CML, given by Eq. (JTJ. 

In general, a chaotically evolving system can be con- 
trolled to a stable periodic orbit by the addition of an 
approp riate constant or periodic external bias 0, El 0, 
IbSl [jjj . As a consequence, one can expect a stable fixed 
point solution for the modified map (jSJ) for appropri- 
ate forcing amplitudes (fci and ^2), with same set of 
parameters for which the original map, x J n+1 — /(x^), 
exhibits chaotic solution. Thus there exists a possibility 
for obtaining periodic solutions of period two for the con- 
stituent map within the CML even though its parameter 
is in the chaotic region of the individual map. 

In the case of coupled logistic lattice with parameters 
mentioned in the previous section, the allowed region 
of forcing amplitudes which lead the decoupled modi- 
fied map © to exhibit period two solution is shown in 
Fig. |21 For L = 8, we have calculated the amplitudes 



of the coupling force [that is, the second term in the 
right hand side of Eq. QJ] of period-2 for the initial con- 
ditions, {x 3 }f^ = {0.1, 0.01, 0.7, 0.2, 0.65, 0.1, 0.15, 
0.001}, and these are shown in Table [D] Now one can 
easily check that these amplitudes of the coupling force 
fall in the region specified by the phase diagram shown 
in Fig. [21 Also, the shift invariant property of the CLL 




FIG. 2: The phase diagram in the ki — plane for the mod- 
ified map © with f(x) = /j,x(1 — x), fi — 3.5732, showing the 
regions which correspond to the periodic solution of period 2 
(dark region). 

ensures that there is no temporal variation if we shift 
the initial conditions of each lattice point in the CML 
to its neighbour spatially. In this case, the wave pat- 
tern will also make only a corresponding shift. We have 
made similar investigations for higher periodic standing 
waves which lead to same type of conclusions based on 
the appropriate periodic nature of the coupling force. 

TABLE II: Amplitudes of the coupling force ki and k% in 
the coupled map lattices for L = 8 with period-2 for initial 
conditions, {x 3 }f~o = {0.1, 0.01, 0.7, 0.2, 0.65, 0.1, 0.15, 
0.001}. 
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So, if it is possible to control the coupling force to fall in 
the region which corresponds to a periodic solution, then 
one can obtain standing wave patterns irrespective of the 
size of the lattices. This is in fact possible by choosing ap- 
propriate initial conditions to each lattice point and this 
explains the occurrence of standing waves (asynchronous) 
as shown in Figs, da) and^b) in the CLL well below 
the critical lattice size (L c ). Similar explanation holds 
good for Fig. ^c) . The same principle is involved in the 
occurrence of standing waves even for L > L c [f| . 
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V. ANALYTICAL EXPRESSION FOR 
STANDING WAVE PATTERNS 

From a careful numerical analysis, we have observed 
that the nodes and antinodes of standing waves are 
formed at or very close to the UPOs of the isolated logis- 
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values from 1 to p and 1 to 2p, respectively and p is 
the period of UPO. Also in the above equation (JJJ, u k 's 
represent the values of the UPOs of the isolated logistic 
map at the node of the fc-th standing wave, and A™ ax and 
A™" are the absolute values of the differences between 
UPOs at the node and UPOs at high and low amplitudes 
at the antinodes of the A:th standing waves, respectively. 
The wave patterns with one and two number of standing 
waves obtained using Eq. J7J) for L = 8 are shown in 
Fig. |3 whereas for the same lattice size, the numerically 
obtained pattern which have been discussed in Sec. Illll is 
shown in Figs. Ufa) and^Jb). One observes that these 
two figures coincide very closely. 



VI. SUMMARY AND CONCLUSION 



FIG. 3: Standing wave patterns obtained from Eq. 10 for 
L — 8: (a) standing wave with period-1 UPO at nodes and 
period-2 UPO at antinodes and (b) period-2 UPO at nodes 
and period-4 UPO at antinodes. 

tic map. Keeping this in mind, we propose an expression 
for the standing wave pattern [2(J of the form 

i a ■ ( mir(j — S)\ , ,. 
x J = u k + A k sin I I cos (m) , (7) 

where, 



In this report, we have pointed out that coupled map 
lattices with diffusicc coupling exhibit multiple stable 
states for the same set of parameters with respect to 
the initial conditions. It has also been shown that by 
choosing appropriate initial conditions, one can obtain 
different standing wave type patterns for such coupled 
map lattices even for lattice size much less than the crit- 
ical system size L c , where one normally would expect 
synchronized chaos. In addition, we have proposed the 
mechanism behind the occurrence of such standing wave 
patterns. 



1 max 

A k J 
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A k J 



if sin (miU—^l \ cos (wi) > 



cos (ni) < 



if sin 
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where the discrete index j — 0,1,2, ■ ■ ■ , L — 1 corresponds 
to the lattice site, m denotes the mode of the waves, and 
k and i represent the nodes and antinodes of different 
standing waves present in the pattern, which can take 
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